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ABSTRACT 

In this paper we present results of a short-term optical monitoring of 13 blazars. 
The objects were monitored mostly in the R-band for a total of ^160 hours between 
2006 and 2011. We study the nature of the short-term variations and show that most 
of them could be described as slow, smooth, and (almost) linear changes of up to ~0.1 
mag/hour, but many objects show no short-term variations at all. In fact, we found 
only ~ 2 per cent chance to observe variability of more than 0.1 mag/hour for the 
sample we observed. Hints for quasi-periodic oscillations at very low amplitude levels 
are also found for some objects. We briefly discuss some of the possible mechanisms 
to generate the intra-night variability and the quasi-periodic oscillations. 
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1 INTRODUCTION 

Nowadays it is widely accepted that a large variety of obser- 
vational blazar phenomena, including the powerful emission, 
the spectral energy distribution (SED) covering practically 
all energy bands, the exceptional variability, the strong po- 
larization, etc. can be attributed to plasma processes in a 
relativistic jet. According to the commonly invoked scheme, 
a highly relativistic jet, pointed at a small angle with respect 
to the line of sight (Urry & Padovani, 1995), generates most 
of the observed spectrum via synchrotron and inverse Comp- 
ton processes. Due to the relativistic beaming, the emission 
is highly anisotropic, significantly amplified, and frequency 
shifted with a typical bulk Doppler factor (8) of around 10- 
30, as assessed by the SED modeling or variability properties 
(Hovatta et al., 2009; Wu et al, 2011b; Rani et al., 2011a). 

As the non jet-dominated AGN, generally do not show 
huge variations, at least not on the intra-night time scales, 
it is almost certain that the blazar variability is also asso- 
ciated with the processes in the relativistic jet (see Wagner 
& Witzel, 1995, for a review). These processes can be in- 
trinsic to the jet and can be related to a change of the jet 
power, for instance. The jet power will change as a result of 
rapid evolution of the energy spectrum of the emitting rela- 
tivistic particles due to either energy loss (both synchrotron 
and Compton losses) or energy gain (re-acceleration or fresh 
particle injection) , or both at different scales. The processes 
responsible for the variations can also be related to the jet 
geometry, i.e. a changing jet direction (a curved or swing- 



ing jet; Copal-Krishna & Wiita, 1992) leading to a vari- 
able bulk Doppler factor of the blobs traveling along the 
jet. The last type of processes generates observable vari- 
ability even if no change of the overall jet power occurs. 
And finally, the variability can be due to entirely extrin- 
sic reasons, like microlensing (Gopal-Krishna & Subrama- 
nian, 1991; Paczynski, 1996) from intervening objects (e.g., 
stars in a foreground galaxy) along the line of sight. Thus, 
the emission from blobs moving relativistically down the jet 
might vary due to microlensing at a much faster rate than 
what is normally expected. Last but not least, the jet itself 
is not a stationary object; a number of reasons including 
turbulence, developing/decaying shocks, different instabili- 
ties, changing environment, etc., can be responsible for the 
variability, at least at longer time scales. 

So, it is not surprising that the exact cause of the blazar 
variability is still under debates. Solving this problem, how- 
ever, is of importance not only to astrophysics. Understand- 
ing blazar variability is also one of the keys to understand 
the physical processes in a relativistic jet. Since the jets are 
natural accelerators to energies, exceeding on many orders 
what can be achieved on Earth, blazar studies could also be 
of extreme importance for understanding the fundaments of 
physics. 

This work focuses on the optical variations on the short- 
est time scales (the so called intra-night or intra-day vari- 
ability). Optical studies are very important for blazars, as 
even though their spectrum is very broad, a key ingredient 
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of their SED, namely the synchrotron peak, is often located 
close to these wavelengths. 

Although time costly, the intra-night variability studies 
cannot be replaced by inter-night (long-term) monitoring, 
since the short-term magnitude gradients are often much 
larger than night-to-night gradients (see for instance Fig. 1 
in Romero et al. 2000; Fig. 2 in Cellone et al. 2007; etc.). 
Fortunately, most of the objects can successfully be mon- 
itored with relatively small (< 1-m) class telescopes, thus 
making possible to organize dense, multisite monitoring of 
selected objects (e.g., WEBT/GASP campaigns^]) . 

In this paper we monitor 13 blazars during 58 nights 
between years 2006 and 2011 for a total of about 160 hours. 
Our goal has been to establish the general character of the 
variations on intra-night time scales. Even though there are 
a number of reports for very rapid, magnitude-scale "frame- 
to-frame" variations, our observations indicate mostly slow 
trends or wobbles that can successfully be fitted with a low- 
order polynomial (Sect. 3), and for many objects/nights we 
find no variations at all. 

We continue next with describing our observational 
data. 



2 OBSERVATIONAL DATA 

We employed several instruments equipped with CCD cam- 
eras and standard UBVRI filter sets to perform the photo- 
metric monitoring. These include the 60cm reflector of Be- 
logradchik Observatory, Bulgaria, (coded B60 in Table 1), 
equipped with SBIG ST-8 CCD (as of 2009 replaced with 
FLI PL9000); the 2m RCC reflector (Photometries AT200 
CCD) and the 50/70cm Schmidt camera (SBIG ST-8, later 
replaced with FLI PL16803 CCD) of Rozhen National Ob- 
servatory, Bulgaria, (coded R50/70 and R200 respectively); 
and the 104cm reflector of ARIES, Nainital, India (Wright 
2k CCD), coded as A104. Further details on the instruments 
in use can be found in Gaur et al. (2012b). 

A total of almost 3000 individual frames were obtained 
and analyzed^. The typical exposure time was 120s; how- 
ever, depending on the brightness of the object and the 
telescope size, different exposure times between 60 and 240s 
were used on different occasions. The monitoring was per- 
formed mostly in the R-band with a few exceptions, where 
the I-band or no filter were used (Table 1). The latter was the 
case for two very weak objects, monitored with a relatively 
small telescope; see Bachev et al. (2005) for comments on 
this issue. The majority of the nights were clear and stable 
but on a few occasions we monitored during non favorable 
atmospheric conditions, which resulted in larger and vari- 
able photometric errors (see Sect. 4.1 for discussions). 

After bias, dark current (where appropriate) and flat 
field corrections, the magnitudes were extracted by apply- 
ing standard aperture photometry routines. The aperture 
radius was taken to be typically 2-3 times the seeing. The 
magnitudes of the blazar and a check star were measured 
with respect to a main standard of known magnitude with 
no further corrections made. Both, the main standard and 

1 Details at http://www.to.astro.it/blazars/WEBT 

2 Some of these data may have been also used in active WEBT 
campaigns 



Table 1. Observations 



Object 


UT date 


Telescope 


Duration 


Points 


Filter 


3C 66A 


14.10.2007 


R50/70 


1.1 


18 


R 


(0219+428) 


02.11.2007 


B60 


2.8 


79 


R 




03.11.2007 


B60 


1.9 


53 


R 




08.10.2009 


A104 


1.6 


31 


R 


1ES 0229+200 


09.10.2009 


A104 


5.3 


67 


R 




21.11.2009 


A104 


4.4 


18 


R 


AO 0235+16 


20.02.2007 


R50/70 


0.9 


20 


R 




25.02.2007 


R50/70 


1.0 


20 


R 


4C 47.08 


10.10.2009 


A104 


3.8 


19 


R 


(0300+470) 


22.11.2009 


A104 


1.1 


13 


R 


S5 0716+714 


16.01.2007 


B60 


1.3 


19 


R 




18.01.2007 


B60 


0.9 


19 


R 




24.02.2009 


B60 


6.6 


161 


1 




25.02.2009 


B60 


6.4 


180 


1 


PKS 0735+178 


08.10.2009 


A104 


1.6 


32 


R 




20.12.2009 


A104 


6.9 


115 


R 


PKS 0736+017 


29.02.2008 


B60 


5.5 


77 


NO 




03.03.2008 


B60 


3.5 


40 


NO 




25.01.2011 


B60 


3.3 


44 


NO 




25.04.2011 


B60 


0.8 


21 


NO 


OJ 287 


19.11.2006* 


B60 


0.9 


26 


R 


(0851+202 ; 


14.lz.2U06 


R200 


4.2 


98 


R 




17.12.2006* 


R50/70 


1.9 


30 


R 




20.02.2007* 


R50/70 


1.9 


30 


R 




08.04.2007* 


R200 


1.8 


10 


R 




09.04.2007 


R200 


2.1 


50 


R 




10.04.2007* 


R50/70 


2.9 


45 


R 




29.02.2008* 


R50/70 


2.8 


75 


R, 


3C 279 


26.01.2006 


B60 


0.5 


15 


R 


(1253-055) 


27.01.2006 


B60 


1.1 


30 


R 




28.01.2006 


B60 


0.7 


20 


R 


3C 345 


20.06.2007 


R200 


3.4 


69 


R 


BL Lac 


15.07.2007* 


R200 


0.4 


10 


R 


(2200+420) 


16.07.2007 


R200 


3.1 


47 


R 




14.08.2007 


R200 


1.9 


50 


R 




15.08.2007 


R200 


2.7 


70 


R. 




16.08.2007 


R200 


2.3 


75 


R 




18.08.2007 


T~> trn /td 

H50/70 


3.4 


70 


R 




19.08.2007 


R50/70 


3.2 


58 


R 




on no n n n t 

20.08.2007 


R50/70 


3.0 


55 


R 




no 1 i ann7 
U2.11.ZUU7 


B60 


0.7 


20 


R 




n""7 1 o onnT* 

07.12.2007 


B60 


2.0 


47 


R 




no n -i onno 

08.01.2008 


B60 


2.6 


61 


R 




09.01.2008* 


B60 


2.3 


58 


R 




11.01.2008 


B60 


1.8 


49 


R 




07.07.2008 


R200 


1.5 


74 


R 




08.10.2009 


A104 


1.0 


78 


R 


3C 454.3 


28.07.2008 


B60 


1.9 


54 


R 


(2251+158) 


29.07.2008 


B60 


2.2 


62 


R 




30.07.2008* 


B60 


2.2 


63 


R 




31.10.2010* 


R50/70 


5.0 


27 


R 




01.11.2010* 


R50/70 


5.8 


27 


R 




02.11.2010* 


R50/70 


5.3 


30 


R 




03.11.2010* 


B60 


6.7 


57 


R 




04.11.2010* 


B60 


3.1 


29 


R 




06.11.2010* 


R50/70 


3.6 


21 


R 


B2 2308+341 


10.08.2010 


B60 


2.8 


40 


NO 



Intra-night variability in blazars 3 



14.10.3007 



03.1 1.3007 




FKS 0736+017 




FKS 0736+017 



PKS 0?:i!-; : I V 



16.3 f-j 
16.4 

16.5 F- 

14.3 3 

14.4 r 

14.5 r 

14.6 - 



Figure 1. Intra-night blazar LCs (3C 66A data). The upper 
panel of each box shows the blazar LC (filled symbols), measured 
with respect to a main standard, and the lower one - the same for 
a check star (open symbols). The photometric errors at 1-rr level 
are indicated as error bars. The best fitting low-level polynomial 
(see the text) is also shown for both - the blazar and the check 
star LCs. The name of each monitored object and the UT date 
when the observation started are shown atop of each box. The 
abscissa shows the UT in hours. Each upper panel has the same 
vertical scale factor as the corresponding lower panel, but the 
exact value may differ from object to object. 
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Figure 3. The same as Fig. 1 (PKS 0736+017 data) 
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Figure 4. The same as Fig. 1 (1ES 0229+200, AO 0235+16, and 
4C 47.08) 
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Figure 2. The same as Fig. 1 (S5 0716+714 data) 



the check star were chosen on the basis of their magnitude, 
color, and spatial proximity to the blazar; however, due to 
the different instrument characteristics (fields of view), no 
formal criteria were applied. 

Table 1 gives a short log of the observations. The 
columns give the object name, the UT date when the obser- 
vation started, the instrument, the monitoring duration in 
hours, the number of the observational points for the corre- 
sponding run, and the filter. Asterisk after the date indicates 
an observation, used for the statistics but not shown in the 
figures. 



3 RESULTS 

Figs. 1-8 show the results of our intra-night variability 
study. The top panel of each box shows the blazar light 
curve (LC) and the bottom one - the check star LCs, both 
measured with respect to the main standard. The object 
name and the UT date are given atop of each box. 

3.1 The nature of the intra-night variations 

As one can see from Figs. 1-8, the blazar LCs can be 
described as smooth fluctuations or trends without any vi- 
olent, " frame-to- frame" changes. Furthermore, no statisti- 
cally significant variations are present at all on many occa- 
sions. Therefore, it would be appropriate to fit the LCs with 
a straight line or - less often (~40% of the cases) - with 
a low-order polynomial. This approach enables us to better 
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Figure 5. The same as Fig. 1 (PKS 0735+178, 3C 279, and 3C 
345) 
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Figure 7. The same as Fig. 1 (BL Lac, larger telescopes data, 
6 out of 7 nights shown) 
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Figure 6. The same as Fig. 1 (OJ 287, 3C 454.3 and B2 
2308+341). For OJ 287 only 2 out of 8 LCs are shown. For 3C 
454.3 only 2 out of 9 LCs, are shown; see also Bachev et al., 2011, 
for the rest of them. 



represent the variations and to eliminate to a large extent 
any spurious effects introduced by the photometric errors. 
Fitting a low order polynomial was already successfully ap- 
plied by Bachev et al. (2011) to search for time delays be- 
tween different optical bands (see also Montagni et al. 2006, 
who fitted straight lines to fit different segments of the LCs). 
The corresponding best fits are also shown in Figs. 1-8 for 
the LCs of both - blazars and check stars. As one sees, the 
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Figure 8. The same as Fig. 1 (BL Lac, smaller telescopes data, 
6 out of 8 LCs shown) 



fits are rather successful in most (but perhaps not all!) cases 
even though some of the runs lasted for more than 5 hours. 
Fitting LCs allows studying the distribution of magnitude 
change rates, dm /At, for all objects individually, as well as 
for the entire sample (Sect. 3.3). 

The order of the polynomial was chosen to be as low 
as possible and a linear slope was used in almost ~ 50% of 
cases. Unfortunately, there is no easy way to use statistics 
(e.g., x 2 ) to determine the degree of the polynomial. The 
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reason is that the real photometric errors are often larger 
(e.g., Bachev et al., 2005, and the references therein) than 
the theoretical ones by a factor that is difficult to determine 
and may well vary from instrument to instrument and/or 
observing conditions. Furthermore, correcting the theoret- 
ical errors by a certain factor does not seem to be a good 
practice (de Diego, 2010). We have to mention, though, that 
the distributions of the magnitude change rates depend very 
little on the exact polynomial order, as long as this order is 
kept much lower than the number of data points. 

Even though we find no rapid variations (say more than 
a few hundreds of a magnitude on a minute time scale) , there 
are sometimes such claims in the literature, some of which 
concern objects from our sample. It is difficult to judge on 
the reality of these results as often there is no independent 
way to verify them. See Cellone et al. (2007) for examples 
and a critical analysis. 

3.2 Individual notes 

3C 66A. During 4 nights of monitoring (Fig. 1), this TeV 
blazar showed only minor (~ 0.05 mag/hour) variations at 
best. This behavior is consistent with the reports by Raiteri 
et al. (1998), Dai et al. (2001), Sagar et al. (2004), Bottcher 
et al. (2009), Rani et al. (2011b), where 3C 66A was observed 
at a similar brightness stage and showed similar short-term 
variability. Although a rapid change was reported on one 
occasion (AB ~ 0.5 mag/hour; Zhang et al., 2004), overall, 
the object does not seem to be highly active most of the 
time, showing only minor, smooth variations. 
1ES 0229+200. We monitored this VHE 7-ray source for 
2 nights for a total of ~ 10 hours (Fig. 4). To the best 
of our knowledge, no results of previous intra-night optical 
monitoring of this object were published, and Kurtanidze 
et al. (2004) found no variations on longer time scales. Our 
monitoring, however, indicates the presence of rapid intra- 
night variability of up to 0.15 mag/hour. As a matter of 
fact, this is one of the few objects, where even a 3th degree 
polynomial did not seem to fit the LC very well. For this 
particular object we had to use two different polynomials to 
fit the entire LC. 

AO 0235+164. We monitored this object as a part of an 
active WEBT campaign (Raiteri et al., 2008a) during 2 pho- 
tometrically poor nights (Fig. 4) and found practically no 
intra-night variations, but a change of AR ~ 0.4 mag was 
seen within ~ 5 days. Other authors (Hagen-Thorn et al., 
2008), however, report significant intra-night variations of 
AR ~ 0.5 mag within several hours. Romero et al. (2000, 
2002) also observed rapid changes, indicating that variations 
of AR ~ 0.1 mag/hour may be typical for this otherwise 
highly active object (see also Sagar et al., 2004; Gupta et 
al., 2008; Rani et al., 2011b). 

4C 47.08. This Fermi LAT source (Abdo et al., 2010) has 
been rarely studied at optical wavelengths. We found a ~ 
0.5 mag variation on a month time scale, but only minor 
variations during the 2 nights of intra-night monitoring. The 
object was rather weak and the photometric uncertainties 
might have contributed to the apparent variations, seen as 
smooth fluctuations (Fig. 4). We note that Star C5 (the 
closest standard, Fiorucci et al., 1998) appears to be variable 
on short time scales. 

S5 0716+714. This blazar has been a target of a number 



of intra-night monitoring campaigns. Practically all stud- 
ies, e.g., Sagar et al. (1999), Nesci et al. (2002), Wu et al. 
(2005), Gu et al. (2006), Montagni et al. (2006), Pollock et 
al. (2007), Gupta et al. (2008), Stalin et al. (2009), Poon et 
al. (2009), Carini ct al. (2011), Rani et al. (2011b), report 
significant variations, which can typically be described as 
smooth fluctuations or a sequence of trends, often abruptly 
changing direction. Some reports even imply the presence 
of very rapid (tenths of a magnitude) changes within a few 
minutes (Fan et al., 2011) and/or quasi-periodic oscillations 
(Gupta et al., 2009; Rani et al., 2010) on ~ 15 min time 
scale. Clearly, S5 0716+714 is one of the most active blazars 
in terms of rapid intra-night variations. Our monitoring dur- 
ing 4 nights showed smooth trends, reaching up to ~ 0.05 
mag/hour but no rapid changes on a minute time scale (Fig. 
2) 

PKS 0735+178. This is another object very extensively 
studied in the optical wavelengths throughout the years. 
Some of the more recent studies by Bai et al. (1999), Sagar 
et al. (2004), Stalin et al. (2004), Gupta et al. (2008), Goyal 
et al. (2009), Rani et al. (2011b) report only low-level intra- 
night variations, if any at all. On the other hand, Zhang et 
al. (2004) find variability of ~ 0.5 mag/hour in one occasion. 
Our observations confirm what most of the researchers re- 
port, i.e. we found practically no variability during 2 nights 
for a total of about 8.5 hours of monitoring (Fig. 5). 
PKS 0736+017. This FSRQ is known to show occasionally 
exceptional variability on intra-night scales. Clements et al. 
(2003), for instance, reported a 1.3 mag outburst within 2 
hours. Their observations, however, as well as our results 
imply that PKS 0736+017 is much more active, when close 
to its maximum (R~13.7 mag), or at least that periods of 
intra-night activity are followed by periods of almost con- 
stant brightness. Sagar et al. (2004) also reported very qui- 
escent behavior when the object was around R~15 mag. 
Our observations (Fig. 3) were during a deep minimum - 
~16 mag (see also Ramirez et al., 2004), yet some low- level, 
smooth variations are seen. 

OJ 287. This object manifested very little if any intra-night 
variability during 8 nights of monitoring, spanning over ~ 15 
months, though night-to-night variations of ~0.1 mag were 
present. In Fig. 6 we only show the results from 2 atmospher- 
ically stable nights, where gradual trends are clearly seen. In 
spite that we found no rapid variations during our monitor- 
ing of this object, other authors report significant changes 
on very short time-scales. Zhang et al. (2007), for instance, 
found a AR ~ 0.35 mag change for a half an hour during one 
night (see also Xie et al., 2001). On the other hand, Carini 
et al. (1992), Gonzalez-Perez et al. (1996), Dultzin-Hacyan 
et al. (1997), Sagar et al. (2004), Gupta et al. (2008), Rani 
et al. (2011b), report only minor variations, if any, which is 
much more consistent with our findings. 
3C 279. The object was monitored in 3 consecutive nights 
in Jan. 2006, being around 15 mag (much fainter than its his- 
torical maximum, R~ 12.5 mag). A clear decreasing trend of 
0.15 mag/hour, probably tracing the end of a small outburst 
(Fig. 5, see also Bottcher et al., 2007, for a long-term LC 
covering this period) was seen in one night, while during the 
other two nights the object was quite stable. Kartaltepe & 
Balonek (2007) reported similar gradual trends (if any) dur- 
ing their monitoring, when the object was around 14 mag. 
Gupta et al. (2008) found practically no variability during 
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their monitoring, but on the other hand, rapid changes were 
reported by Miller et al. (2011) for low state. Clearly, 3C 279 
exhibits different types of intra-night variability behavior. 
3C 345 This superluminal FSRQ has not been observed 
for intra-night variability very often, probably because it 
has been rather weak (~ 16 — 17 mag) in the recent years. 
Mihov et al. (2008) and Howard et al. (2004) found practi- 
cally no intra-night variability. On the other hand, Wu et al. 
(2011a) reported trends of ~ 0.2 mag/hour, while Kidger & 
de Diego (1990) found even sharper decline in one occasion. 
We monitored this object during one night and found no 
variations at all (Fig. 5). 

BL Lac. Wc monitored extensively the blazar archetype for 
a total of almost 35 hours during 15 nights, in 6 of which 
a 2-meter class telescope was used (Figs. 7, 8). Variability 
pictures reveal trends of up to ~0.1 mag/hour or smooth 
fluctuations. The object was mostly around R ~ 14.5 mag, 
meaning a rather low state, and showed little night-to-night 
changes (Raiteri et al., 2009). Other researchers report sim- 
ilar intra-night behavior (Nesci et al. 1998; Speziali & Na- 
tali 1998; Papadakis et al. 2003; Howard et al. 2004), e.g., 
smooth variations of typically ~0.1 mag/hour. On the other 
hand, Zhang et al. (2004) found unusual ~0.3 mag frame- 
to-frame variations. 

3C 454.3 This is another well studied FSRQ. Some of 
the recent intra-night variability studies include Poggiani 
(2006), Gupta et al. (2008), Rani et al. (2011b), Gaur et al. 
(2012a) who found only minor or no variations, but ~ 0.2 
mag/hour variability rates were seen occasionally. Similar 
behavior (see also Fig. 6) was also reported by Bachev et al. 
(2011), who monitored the object during a recent outburst 
to search for time delays between different optical bands on 
intra-night time scales. Their R-band LCs were also used 
for the purposes of this work. Zhai et al. (2011) performed 
a similar search but none of these works found convincing 
evidence for a delay between the optical bands. Raiteri et 
al. (2008b), on the other hand, give clear examples of vi- 
olent intra-night variability episodes of ~ 1 mag/hour per- 
haps tracing short-lasting outbursts during the high activity 
phase of the blazar. Obviously, this object demonstrates dif- 
ferent variability behavior, probably depending on its bright- 
ness state. 

B2 2308+341 Optical monitoring of this object was trig- 
gered by increased 7-ray activity (D'Ammando, 2010). The 
object was very weak, so no filter was used during the 2 
nights of monitoring. Thus, the magnitudes shown in Fig. 6 
are somewhat arbitrary. No huge intra-night or day-to-day 
variations were observed. The object has been rarely studied 
and we found no other works on intra-night variability. 

3.3 Distribution of the variability rates of the 
sample 

The LCs we analyze in this work cover a total of ~ 160 hours 
of monitoring. Although the LC data can not be considered 
complete nor can the sample be considered representative, 
it still may be worth building a histogram of the magni- 
tude change rates (dm /At) for the entire sample. Such a his- 
togram may help studying time asymmetries, for instance. 
Fig. 9 shows the distribution of dm /At values for the blazars 
(upper panel) and the check stars (lower panel) based on 
the corresponding best fits. Normally, one would not expect 
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Figure 9. Magnitude change rates for the blazars and the check 
stars, as found from the fitting polynomials. The numbers in the 
distribution count the 3 min intervals from the LC with the cor- 
responding magnitude change rate. 

the check stars to show variations; however, sometimes slow 
trends became evident after applying the same fitting proce- 
dure as for the blazars, resulting in some non-zero values in 
the dm /At histogram for the check stars. Explanations can 
be searched for in terms of variable atmospheric conditions; 
possible low-level stellar variability; etc. (see Klimek et al., 
2004; Bachev et al., 2005), but this effort goes beyond the 
scope of this paper. 

There are two important results that are worth men- 
tioning. First, it seems that observing significant variabil- 
ity rates is a very rare case. For instance, from Fig. 9 
we find only ~ 2% chance to observe blazar variability 
|dm/dt| > 0.1 mag/hour and ~ 25% chance to observe prac- 
tically no detectable variability, |dm/dt| < 0.005 mag/hour 
(the corresponding percentage for the check stars is ~ 68). 
Secondly, the distribution does not appear to be symmetric. 
Statistical tests to compare the means, medians, standard 
deviations, as well as the KS test, performed separately on 
the "negative" and "positive" parts of the Am /At distribu- 
tion distinguish them at 99.9% level, even though the mean 
of the entire distribution is consistent with zero within the 
errors (as expected). This might be a result of the data in- 
completeness; however, it also might be an artifact of a real 
time asymmetry of the blazar LCs. Since the total monitor- 
ing times were very different for different objects, we cannot 
claim that such asymmetry, if real at all, is a general prop- 
erty of the sample or, at even lesser extent, of the entire 
blazar population. 

3.4 Possible quasi-periodic oscillations 

Employing a 2-m class telescope to monitor some of the ob- 
jects (BL Lac in particular) ensured photometry of very high 
accuracy (e.g., 0.002 - 0.003 mag) during some stable nights, 
which on the other hand allowed studying possible (quasi- 
periodic) micro-oscillations. Such have been previously re- 
ported; e.g., Clements et. al. (2003) for PKS 0736+017; Wag- 
ner et al. (1996), Gupta et al. (2009), Rani et al (2010) 
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Figure 11. The influence of the atmospheric transparency on the 
differential light curves of the blazars - S5 0716+714 and PKS 
0736+017 monitored during two unstable nights. The upper boxes 
show the blazars' LCs and the lower ones - the signal from the 
main standard, transformed into arbitrary magnitudes. One sees 
that the atmospheric stability does not have any visible effect on 
the differential LCs. 



Figure 10. BL Lac light curves after subtracting the leading 
polynomial (Fig. 7). Moving average filter is applied to the data 
to reduce the scatter. Micro-oscillations are clearly seen in some 
nights (e.g., 14.08.2007). The upper panels show the sine func- 
tion best fits, and the lower ones — the residuals. For the night 
of 15.08.2007, an exponentionally decaying sine amplitude works 
best. 

for S5 0716+714; etc. For BL Lac, some hints for quasi- 
periodicity were presented by Speziali & Natali (1998, their 
Fig. 3), though the authors do not claim any statistical sig- 
nificance. It is clear from our BL Lac data (the most exten- 
sively studied object in our sample) that a periodic signal is 
not the dominating feature in the LCs. However, after sub- 
tracting the leading polynomial, one may attempt to search 
for micro-oscillations in the residuals. We found some clues 
for the presence of such micro-oscillations with an amplitude 
of ^ 0.01 mag and a typical period T ~1 hour for many of 
the nights (Fig. 10). 

Interestingly, the oscillation pattern does not appear to 
depend much on the order of the leading polynomial (within 
some limits, of course), used to fit the LC. If further con- 
firmed, micro-oscillations with a period of an hour or so may 
have significant implications for our understanding of how 
matter accelerates along the jet. We note on the other hand, 
that such a period is close to the shortest possible Keplenan 
time (perhaps even below the innermost stable orbit, Sect. 
4.3). 



4 DISCUSSION 

4.1 Differential light curves and atmospheric 
stability 

Since the colors of the blazars and the stars used for the 
differential photometry may sometimes differ significantly, 
one can raise the question about the influence of the atmo- 
spheric stability and/or air mass change on the differential 
LCs. In order to check for such effects we consider closely 
cases when the blazar showed significant variations during 
the monitoring in unstable nights. Fig. 11 shows two such 
cases. The upper boxes show the blazars' LCs and the lower 
ones - the signal from the main standard, transformed into 
arbitrary magnitudes. One sees that while the main stan- 
dard's magnitude varied significantly, both because of the 



changing air mass and the unstable atmosphere, the differ- 
ential blazar LCs do not seem to be influenced at all, which 
leads to the conclusion that the observed trends are un- 
likely to be attributed to the changing atmospheric trans- 
parency. The Pearson correlation coefficients between the 
two datasets are -0.36 for S5 0716+714 and -0.66 for PKS 
0736+017, implying a rather insignificant correlation. Spear- 
man rank correlation coefficients are even less significant - 
—0.27 and —0.41, respectively. As a matter of fact, a small 
anticorrelation between the datasets is actually expected, as 
the photometric errors of the main standard influence both 
datasets. 

4.2 Implications for the variability models 

4-2.1 Time asymmetry of the LCs 

Studying the magnitude gradients' distribution may help 
to disentangle among different variability scenarios. For in- 
stance, a dissipative process (i.e. an explosive event, etc.), 
being by nature not time reversible, will likely produce a 
time-asymmetric LC (perhaps with sharp raises and slow 
declines), which will result in different positive and negative 
magnitude change rate distributions. On the other hand, 
a process described by a conservative dynamical system, 
whose Hamiltonian is not directly depending on timqj (like 
microlensing, orbital motion, etc.) will likely produce a time- 
symmetric LCs, leading to a symmetric positive/negative 
gradient's distribution. Unfortunately, our data (Fig. 9) can- 
not be considered statistically significant for a number of 
reasons (Sect 3.3) for a firm conclusion. 

4-2.2 Evolution of the electrons' energy density 
distribution 

A natural cause of the brightness changes can be the evolu- 
tion of the relativistic particles' energy density distribution, 
71(7), which consequently leads to a variable synchrotron 
emission. While the energy loss rates are known (both syn- 
chrotron and Compton losses are oc j 2 , 7 is the electron's 
Lorentz factor; Rybicki & Lightman, 1979), the energy gain 
rates (the acceleration/injection mechanisms) are generally 

3 See however for more details, e.g., Lamb & Roberts, 1998, and 
the references therein 
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Figure 12. A sketch showing the change of the jet direction and 
the corresponding angles (see the text). 

unknown. This leads to the impossibility to solve the cor- 
responding Fokker-Planck equation (Rybicki & Lightman, 
1979) without a priori imposed assumptions. Electron en- 
ergy density evolution will probably lead to time delays be- 
tween different wavebands (Bachev et al., 2011, and the ref- 
erences therein). Knowing that our LCs are monochromatic 
and taking into account the above said, places testing this 
variability mechanism beyond the scope of this paper. 

4-2.3 Swinging jet 

The idea behind this variability mechanism is that the path 
of the relativistically moving blobs along the jet can devi- 
ate slightly from a straight line, thus leading to a variable 
Doppler factor and respectively to variable emission from 
that particular blob (Fig. 12 ilustrates the situation, see 
also Gopal-Krishna & Wiita, 1992). From the definition of 

the Doppler factor 8 = ^ — | ^ = — — t T where 

r(i— p cos 8) r(i — pi^o cos (fij ' 

/io = cos#o and 9o is the "impact parameter", i.e. the small- 
est angle between the jet direction (direction of motion of 
the blob) and the line of sight. Clearly, the maximum am- 
plification occurs when 9o = 0. 

If take into account that all angles are small (we con- 
sider only cases where the magnification is significant), it 
can be shown that S = , , A ..> , where A is a constant and 
t is a characteristic timescale. 

The magnitude change associated (if this blob is dom- 
inating the blazar emission at some particular moment) 
can be calculated taking into account that the flux F v = 
i 7 ^/ <5 3+a , where a is the spectral index, i.e., F v oc v a . Thus, 
dm/dt ~ dln<5/di (see also Nesci et al., 2002; Montagni et 
al., 2006). Knowing dm/di one can integrate to obtain the 
magnitude of the blob emission as a function of time. 

As expected, m(t) is in general a complex peak func- 
tion, but for practical purposes it can be approximated with 
a simpler one. For t >> r, the peak appears quite sharp 
and the best working simple function here is the exponent. 
However, the magnitude change becomes significant; e.g., 
~ 10 mag (see also Gopal-Krishna & Wiita, 1992). Since 
such large amplitudes are never observed in the blazar LCs, 
one has to assume that a flare from a single blob, due to 
a variable Doppler factor, may never completely dominate 



the blazar emission, except perhaps for t ^ r. Thus, a cor- 
rect approach would be to invoke an additional constant flux 
component of magnitude mi , which can be due to the host 
galaxy emission and/or to an ensemble of independent flar- 
ing events, each at different evolutionary state (but yet, the 
most of them quite far from the maximum). 

Introducing the constant component changes signifi- 
cantly the shape of the flares; this time the best working fit 
appears to be a Gaussian for a fairly broad range of ratios 
between m and mi. In other words, unless we accept that 
a flare can be so powerful to increase the total blazar emis- 
sion by many orders (say by more than 10 mag), its overall 
shape has to resemble a Gaussian within the framework of 
this model. 

4-2.4 Microlensing 

We can use a similar scheme to compute the LC of an emit- 
ting blob, which emission is magnified as a result of mi- 
crolensing from a star (a point-mass lens) along the line of 
sight (Paczynski, 1996). The flux magnification for gravita- 
tional lensing is M = — °, +2 , where 8 — 9/9e, where 8e 

is the Einstein radius, (see for instance Narayan & Bartel- 
mann, 1996, for a review). 

Under proper assumptions one can calculate M(t) and 
respectively - the light curve if the emission is dominated by 
a single, microlensed blob. In this case the best working fit 
appears to be the Lorentzian, gradually being transformed 
into an exponent for strong microlensing events, see also 
Paczynski (1996). Interestingly, adding a constant flux as 
before does not seem to affect significantly the shape of the 
flare (within some limits). 

4-2.5 Comparison with observations 

A dense LC covering a relatively long time interval may 
allow studying the shapes of the individual flares and thus 
help to constrain the models. As we saw, the latest two mod- 
els we considered in detail, predict in general different flare 
shapes, which is especially true for the stronger flares. 

On the other hand, a number of authors have success- 
fully modeled the blazar flares with exponents (Valtaoja et 
al., 1999; Bottcher & Principe, 2009; Hovatta et al., 2009; 
Chatterjee et al., 2011), which resembles the case of a strong 
microlensing evenlQ. Most of these authors, however, used 
different rising and decaying times for their exponential fits, 
which cannot be easily explained in terms of microlens- 
ing (see also Reinthal et al., 2012, for statistics on time- 
asymmetry for many objects) . If not attributed to a possible 
degeneracy of the fitting procedure when applied to multi- 
ple flares, such a time asymmetry might be a result from 
the evolution of the electron energy density distribution, i.e. 
from different characteristic times of the acceleration and 
energy loss mechanisms. Having said that, we stress again 
that to calculate the emission one needs to solve the Fokker- 
Planck equation. For example, even if the energy loss is the 
only mechanism to modify the energy density distribution 

4 One is to bear in mind, however, that many of these authors 
used flux units, not magnitudes, when fitting the flares 
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of the relativistic particles, this would not necessarily mean 
a decreasing flux at certain wavelengths and vice vers4f|. 
Actually, the asymmetric flares can also be a result from 
jet inhomogeneity with a combination of jet curvature (e.g., 
Villata et al. 2009). 

Unfortunately, during our monitoring we could not 
trace individual flares, so it is difficult to distinguish which 
one (if any) of the models is mostly responsible to produce 
the variability (the flares) or even which one plays more im- 
portant role. The situation becomes even more complicated 
if one considers a combination of different flares with dif- 
ferent phases and amplitudes to produce the observed vari- 
ability picture. In any case, the presence of very sharp peaks 
seems to favor microlensing. It is also possible that different 
mechanisms govern the flux changes at different time scales. 

4.3 Radiation from the black hole vicinity? 

Possible quasi-periodic oscillations with a period of an hour 
or so (the BL Lac case, Fig. 10), if associated with Keple- 
rian motion around the supermassive black hole, may im- 
pose certain constraints on the black hole characteristics - 
the mass (Mbh) and the spin parameter (a). The Keplerian 
period (in hours) at a distance r is tn — (r 3//2 +a)Ms, where 
Ms is the central mass in 10 8 Solar masses (e.g., Shapiro 
& Teulkolsky 1983), and r is in units of M (geometrical 
units used, G = c = 1). Since the horizon is determined by 
rn = 1 + (1 - a 2 ) 1/2 , one gets (taking that t K ^ 1) M s < 0.5 
for any value of the black hole spin (0 ^ a ^ 1). The mass 
should be even smaller if one accepts a much more physical 
assumption, i.e. that the oscillation is due to some instabili- 
ties at the innermost stable circular orbit (nsco), instead of 
at the horizon. The last distance depends in a complex way 
on the spin and the mass (e.g. Shapiro & Teulkolsky, 1983). 
A simple useful fit for practical purposes that works well 
everywhere, except for a ~ 1, and combines tK (in hours) at 
risco with the black hole parameters is the following expres- 
sion: M 8 t K /V330.7 e~ a - 115.2. Therefore, M BH should 
be between 7 10 6 M Q for a = and 4 10 7 M Q for a = 0.998 
(Fig. 13) if the oscillations are due to Keplerian motion at 
nsco- Such a mass is at least an order of magnitude lower 
than what the host galaxy magnitude suggests for BL Lac 
(Wu et al. 2002; Falomo et al. 2003), but is similar to other 
estimates, based on the variability characteristics (Xie et al. 
2002). It is clear, that since the link between these possible 
quasi-periodic oscillations and the orbital motion is far from 
proven, no firm conclusion is possible at this stage. Studying 
such oscillations, however, if furhter confirmed with high- 
accuracy observations, can become one of the very few ways 
to probe the blazars' innermost regions. 



5 CONCLUSIONS 

We publish optical intra-night LCs of 13 blazars (BL Lacs 
and FSRQs), covering a total of ~ 160 hours of monitor- 
ing time. In the majority of cases we detected no variability 
at all, and when the variations were found, they resembled 

5 If due to energy loss, the synchrotron flux will drop at all wave- 
lengths only if a > 2, where 71(7) oc "f~ a 
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Figure 13. The radial distances for an orbital period of 1 hour as 
functions of black hole mass and spin parameter. Thick continuus 
lines indicate the innermost stable orbit and the horizon, dashed 
lines — the 1 hour period distances for Mg = 0.07, 0.1, 0.25 and 
0.4, from top to bottom, respectively. Note that all the solutions 
passing below the horizon line are unphysical there. 

mostly linear trends or smooth fluctuations. In any case, we 
observed no erratic, "frame-to- frame" changes, which are 
sometimes reported. The smoothness of the LCs allowed fit- 
ting a low-order polynomial, which helped to better reveal 
the magnitude trends. 

The distribution of the magnitude change rates ap- 
pears to be time-asymmetric, which, if real, would mostly 
favor models, based on electron energy density evolution 
rather than "time-symmetric" models like jet swinging or 
microlensing. 

We found some indications for the presence of quasi- 
periodic micro-oscillations with a period of about an hour in 
the high accuracy LCs of BL Lac after subtracting a leading 
polynomial. 

Unfortunately, at this stage, we were unable to establish 
which one of the existing models is mostly responsible for 
the blazar variability. It is possible that more than one of 
them can play some role for the different objects and/or 
activity episodes of the same object. Clearly, further high 
accuracy intra-night variability studies are needed for a firm 
conclusion. 
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